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Abstract. During the past few years, first observations of starburst galaxies at > GeV ener- 
gies could be made with the Fermi Gamma-ray Space Telescope (GeV range) and Imaging 
Air Cherenkov Telescopes (TeV range). The two nearest starbursts, M82 and NGC253 were 
detected, and most recently, the detection of two starburst-Seyfert composites (NGC1068 
and NGC4945) were reported. The emission for the two starbursts is best explained by 
hadronic interactions, and thus providing a first, unique opportunity to study the role of 
cosmic rays in galaxies. In this paper, the role of cosmic rays for the non-thermal com- 
ponent of galaxies is reviewed by discussing the entire non-thermal frequency range from 
radio emission to TeV energies. In particular, the interpretation of radio emission arising 
from electron synchrotron radiation is predicted to be correlated to TeV emission coming 
from interactions of accelerated hadrons. This is observed for the few objects known at TeV 
energies, but the correlation needs to be established with significantly higher statistics. An 
outlook of the possibility of tracing cosmic rays with molecular ions is given. 
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1. Introduction 

Galaxies with a high star formation rate pro- 
vide the opportunity to study the influence of 
massive stars on the large-scale behavior of 
galaxies in detail. In particular, the role of gas 
heating by star formation can be studied by 
the observation of far-infrared emission and 
the role of supernova remnants can be inves- 
tigated by looking at non-thermal radio emis- 
sion from electron synchrotron radiation. Most 
recently, first gamma-ray detections using the 
Fermi Gamma-ray Space Telescope (FGST) 
and Imaging Air Cherenkov Telescopes like 
H.E.S.S. and VERITAS were made. The near- 
est starburst galaxies M82 and NGC253 were 
detected above GeV energies, tracing the inter- 
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action of hadronic cosmic rays with the gas in 
the galaxy. 

In this paper, the multifrequency spec- 
trum of starburst galaxies will be reviewed. 
In Section |2] the correlation between the far 
infrared and radio emission is discussed. In 
Section [3] hadronic interactions at energies 
E > GeV are discussed in the context of neu- 
tral secondary production. In Section [H cos- 
mic ray-induced ionization is mentioned as a 
future method of tracing hadronic cosmic rays 
in galactic environments. In particular, a su- 
pernova remnant interacting with a molecular 
cloud provides an optimal environment for the 
production of Hj, at a level which should be 
detectable in the future with instruments like 
Herschel and ALMA. 
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2. The Radio-FIR correlation 

The correlation between the total radio and far- 
infrared emission in starburst galaxies is ex- 
perimentally well established. As an example, 
a sample of local starburst galaxies with red- 
shifts below z - 0.03 is shown in Figure [T] 
Far-infrared emission arises when newly born, 
massive stars heat the surrounding gas. Radio 
emission, on the other hand, is directly cor- 
related to the death of massive stars: elec- 
trons are accelerated at supernova remnant 
shock fronts and lose their energy via syn- 
chrotron radiation at radio wavelengths. Thus, 
from this very general statement, a correla- 
tion between radio and far-infrared emission 
is expected. A detailed modeling of the cor- 
relation, however, remains difficult when try- 
ing to explain all observations. While a calori- 
metric model, in which all electrons lose their 
entire energy to synchrotron radiation seems 
feasible i n the t heoret ical framework as pre- 
sented bv lVolk I (Il989l) . observations of the ra- 
dio spectral index in starburst galaxies show 
relatively flat spectra: On average, radio spec- 
tra in the GHz range behave as ~ v '^^. The 
synchrotron spectral index aradio and primary 
particle s pectral index gp are direc tly con- 
nected as (iRvbicki & Li ghtmanll 19791) 



2 ■ ffradio + 1 ■ 



(1) 



Thus, a synchrotron spectral index of cKiadio ~ 
0.7 reveals a primary electron population with 
a spectral index of Up ~ 2.6. In the calorimet- 
ric model, however, higher energy electrons 
should lose all their energy, resulting in spec- 
tra as steep as E^^. An alternative model, based 
on a cosmic ray driven wind model for star- 
burst galaxies, is discussed by Bec ker et al. I 
(l2009h . Here, the non-thermal radio flux, in- 
duced by electron synchrotron losses, turns out 
to be only weakly dependent on the magnetic 
field of the starburst {Fy oc B"'^), even without 
the assumption of a calorimeter. The question 
remains whether or not electrons can partly es- 
cape the starburst galaxies and more detailed 
calculations have to be performed in order to 
match all observational features. 

An additional complication in the interpre- 
tation of electromagnetic radiation from star- 



burst galaxies is the apparent co-existence be- 
tween starburst galaxies and active Seyfert 
cores. Figure \T\ shows those galaxies that re- 
veal an active core in a starburst galaxy as filled 
triangles. Thus, when the emission cannot be 
resolved spatially, ambiguities can arise con- 
cerning the electromagnetic contribution from 
the central activity and the starburst part of the 
galaxy. 

3. Hadronic interactions atE> GeV 

A significant part of the total energy budget 
of a galaxy goes into the acceleration of cos- 
mic rays. It is known from the observation of 
the charged cosmic ray flux at Earth, that the 
total energy carried by Galactic cosmic rays 
with energies above 100 GeV corresponds to 
a total luminosity of ~ lO"*^ erg/s. The shock 
fronts of supernova remnants are one of the 
primary candidates for the accele ration of cos- 
mic r a ys up to 10'^ eV or h igher (St anev et al. I 



mic r a ys up to lU ev ornigner ( at anev et al. I 
Il993t iBiermann et al. Il2009. .2010aiibl) . When 
interacting with the local gas of the consid- 
ered galaxy, a significant amount of the pri- 
mary cosmic rays' energy is put into the pro- 
duction of high-energy photons and neutrinos. 
While neutrino telescopes have not yet reached 
the sensitivity to detect diffuse emission from 
the Milky Way, gamma-ray emission is de- 
tected at a level of lO^** erg/s. Thus, the trans- 
port of cosmic rays in a galaxy is expected to 
play a significant role considering the dynam- 
ics of the galaxy itself. The difficulty in pin- 
pointing the sources of cosmic rays themselves 
lies in the complexity of the transport of the 
charged particles through the galactic magnetic 
field. For the Milky Way, the observed flux of 
cosmic rays contains information about the to- 
tal energy budget and the spectral behavior af- 
ter transport. However, no direct information 
on the direction of the sources of cosmic rays 
or about the primary spectra at injection can be 
deduced from observations. Concerning star- 
burst galaxies, compared to the high cosmic 
ray flux from the Milky Way, any possible sig- 
nal in charged cosmic rays would be negligi- 
ble. So, in general, in order to understand the 
role of cosmic rays in a galaxy, the search for 
neutral secondaries like photons and neutrinos 
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Fig. 1. FIR-radio correlation for a sample of local starburst galaxies at z < 0.03, showing a linear correlation 
between the two wavelengths. Luminosities are given in units of erg/s. The sample contains a significant 
fraction of sources which, apart from the starburst behavior of the galaxy, reveal Seyfert-like properties of 
the galaxy core. While those galaxies with no Seyfert core identification are shown as open triangles, the 
supposed Seyfert- starburst composites are shown as filled triangles. 



from hadronic interactions is one of the most 
interesting approaches. 

High-energy photons and neutrinos are 
produced in hadronic interactions, in environ- 
ments with a high flux of high-energy cos- 
mic rays JpjEp) and a large target density hh 
dfiecker et al. II2009I) : 

PP^N{7T^/7T-/7T°)+X. (2) 

Here, N(n^ jn^ jn^'') denotes that multiple pi- 
ons A^ are created in the process. Charged pi- 
ons contribute to the production of high-energy 
neutrinos. 



n' -^ H"v^^ e' Vf v^ v^ . 



(3) 
(4) 



Neutral pions produce high-energy photons. 

Assuming that the dominant part of 
charged cosmic rays below 10'^ eV is accel- 
erated in supernova remnants, the target den- 
sity directly at the acceleration site is im- 
portant as a local effect. At high hydrogen 
densities, cosmic rays interact directly at the 
source and the neutral interaction products re- 
veal the spectral shape at injection. For sources 
that have proton-proton optically thin environ- 
ments, most protons will escape the acceler- 
ation region and change their energy spec- 
trum due to transport effects. Then, interactions 
with the diffuse gas in the galaxy serve as a 



tracer of the cosmic ray spectrum after trans- 
port. Measurements concerning the transport 
of cosmic rays through the MiUcy Way's mag- 
netic field indicate that diffusion steepens the 
injection spectrum by a factor of ~ E^^-^ — 
E-°-^ (Gupta & Webber 1989). With the ob- 
served spectrum being close to E^^-"^, the ex- 
pected cosmic ray energy spectrum at injec- 
tion is expected to be close to E"^' - £'"^'*. 
At gamma-ray energies above ~ 300 MeV, the 
photon spectrum roughly follows the charged 
cosmic ray spectrum in the case of hadronic 
interactions. Thus, the observation of the spec- 
tral behavior of gamma-rays from a galaxy can 
help to identify if the interactions on average 
take place close to the remnant, at injection, or 
rather in the interstellar medium, after trans- 
port. While observations of our own Galaxy 
indicate that propagation steepens the cosmic 
ray spectrum before interactions take place, the 
two starburst galaxies with observed gamma- 
ray spectra, M82 and NGC253, reveal a quite 
flat spectral behavior of around E^^-^. Thus, it 
is expected that a large fraction of the cosmic 
ray interactions happen in the vicinity of the 
sources. 
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3.1. Gamma-ray observations of 
starburst galaxies 

As of today, there are six objects observed at 
gamma-ray energies, where the emission ap- 
pears to arise from hadronic interactions: 

- M82 (Starburst galaxy) 

- NGC253 (Starburst galaxy) 

- The Milky Way 

- Star-forming region 30 Doradus (Large 
Magellanic Cloud) 

- NGC1068 (Starburst-Seyfert composite) 

- NGC4945 (Starburst-Seyfert composite) 

Ahhough the statistics with only six objects re- 
vealing hadronic interactions is still low, it is 
worth to have a first look at what kind of cor- 
relations to expect and to investigate if those 
show up in this still very small sample. It is 
clear that there are many caveats when looking 
for correlations in this small sample: the only 
true starburst galaxies are M82 and NGC253. 
The star-forming region 30 Doradus is only a 
part of the dwarf galaxy. The Milky Way is 
a regular galaxy with contributions to the FIR 
emission from stars that do not explode as su- 
pernovae. Thus, the Milky Way is not expected 
to lie on the FIR-radio correlation curve. The 
two Starburst-Seyfert composites could have 
contributions from the active core in all wave- 
length ranges. Thus, future observations of a 
larger sample of pure starburst galaxies will 
have to confirm or reject the results presented 
below. 

3.2. Interpretation of measured 
gamma-ray spectra as the product 
of hadronic interactions 

From Fermi- and lACT-observations, the 
gamma-ray flux from M82, NGC253, the 
Milky Way and the star-formin g region 30 
Dorad us in the LMC are known (lAbdo et al. I 
l2010al) . The measurements give a unique op- 



portunity to draw conclusions about the pri- 
mary cosmic ray flux. 

First of all, the gamma-ray spectra for 
M82 and NGC253 appear to be very flat, i.e. 
E^^-^ and flatter In the observed energy region 
(> GeV), the hadronic gamma-ray flux follows 
the spectral behavior of the primary cosmic 
rays. 

Since the gamma-ray spectrum reproduces 
the same spectral shape as the interacting 
primary particles, most of the interactions 
must happen close to the acceleration re- 
gion: while the injection spectru m is expected 
to be close to £~^° - E^ ^-^ ('Stanev et al. I 
[1993 ; Bie rmann et al. I l2009l l2010a b). trans- 
port through the galactic magnetic field steep- 
ens the charged primaries' spectrum to E^^-'' . 
Interaction after transport into the ISM would 
therefore lead to a steep gamma-ray spectrum 
of close to E^^-'', 



dNy 



Jy 



dEy dt dAEarth 



(5) 



in units of GeV" s" cm" . 

The number of pion-decay induced 
gamma-rays of a single SNR in a starburst 

Sier unit tim e, vo lume and energy is given 
Kelner et al. Il2006h . 

% (Ey) = riH ■ 



Je, 



'Ey \ dE„ 

o-mei(£p) ■ j„(E„) ■Fy\-r'Ep\-r- -(6) 



Here, uh is the density of the ambient medium 
and (TiaexiEp) is the cross section of inelas- 
tic proton-proton interaction^. The function 
Fy(x, Ep) implies the number of photons in the 
energy interval (x, x + dx) per collision and 
is a dimensionless probability density distribu- 
tion function. The flux of relativistic protons at 
the source is given as 



JpiE,,) ^ 



■<i>{E,). 



(V) 



Here, jp is given per energy, time and area 
interval. The spectral shape of the interacting 
proton flux is contained in the function 0(£'/>) 
and can be expected to follow a power-law. 



' iKelner et al. I ( l2006h use the cosmic ray density, 
while here the cosmic ray flux is used. 
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the speed of light. Solving the equation for Op 

gives; 



Wn 



(13) 



Thus, the normalization of the photon spec- 
trum can be written as 



(J) ^ £-p^ while Op is defined as the normal- 
ization factor, i.e. how many protons there are 
per energy, time and area interval. This normal- 
ization can be done at an arbitrary energy, we ^ _ 
chose Ep = 1 GeV to normalize the spectra. V ■ f" 'i'(Ep) Ep .dEp 

Using other values yields the same result when 
applying the same units consistently through- 
out the calculation. 

To account for the total gamma-ray flux as 
observed at Earth, in units of per time, area and 
energy interval. Equation Q needs to be multi- 
plied by the volume of the interaction region V 
of a single SNR and the number of SNRs A^snr, 
which scales directly with the SN rate /?sn of 
the galaxy. Further, assuming isotropic emis- 
sion, the detected flux scales with 1/(4 Trii^). 
The observed flux is then directly connected to 
the produced density as 



Wp-c-riH ■ NsNR 
4nd^ ■ J^ <:>(Ep) Ep dEp 



(14) 



®le 



cD^(£^) • V ■ A'sNR 



[Ancfy 



(8) 
(9) 



r°° lEy \ dEp 

■ I cri„ei(£p) • jp(Ep) • ^r — • Ep\ — ■ 

Assuming a proton flux at the source follow- 
ing a power-law with normalization Up and a 
spectral shape (t>{Ep) yields 



Given the distance to the source and determin- 
ing the average shape of the primary particle 
spectrum from the shape of the gamma-ray 
spectrum, the product of the total cosmic ray 
energy budget, the target density and the num- 
ber of SNRs, Wp ■ riH ■ A^SNR determines the 
gamma-ray flux normalization. Alternatively, 
the cosmic ray energy density pcR - Wp/V can 
be used, and in this case, it is the product with 
the total interacting mass. 



PCR • M ■ A^SNR = PCR -riH -V ■ NsNR 



(15) 



I'le 



(10) 



j o-,,,,(Ep) ■<:>„(£„)■ fJ^, E,y^ 



with 



Uy = 



Op - riH-y ■ NsNR 



(11) 



This factor needs to be fixed to a certain value 
to match the observations, and conclusions 
about the parameters on the right-hand side can 
be drawn. 

The proton spectral normalization, ap can 
be calculated using the conservation of energy: 



I ip Ep 



dE„ 



ap- \ <b{Ep)E„dE„ 



w„ 



V 



(12) 



where Wp is the total proton energy budget 
of protons with a minimum energy of E^i„, 
assuming the protons travel approximately at 



is the determining factor Thus, assuming a 
constant cosmic ray energy density and a di- 
rect scaling of the number of SNRs with the 
supernova rate /?sn would result in a direct pro- 
portionality between the observed gamma-ray 
emission and the total interacting gas mass. 

First Fermi results seem to suggest that 
a correlation between the gamma-ray lumi- 
nosity and the product of the superno va rate 
and the total gas mass is present (Abdo et al.l 
l2010ah . However, the calculation has several 
caveats: For once, it is not likely that the cos- 
mic ray energy density is constant for all star- 
burst galaxies. In addition, it is assumed here 
that all SNRs have the same cosmic ray energy 
spectrum, although it is known that the spec- 
trum actually strongly depends on their age and 
the local environm ent of the sources, see e.g. 
iBlasi et aO (l2005l) and references therein. One 
might consider to rather search for a correla- 
tion between the product of the SN rate and 
the gas density with the total gamma-ray flux. 
However, even in this case, the above calcu- 
lations rely on the fact that the energy budget 
for each SNR is the same. However, it is ex- 
pected that it depends on the total mass of the 
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stars and the energy put into cosmic rays might 
not even be a constant for the same progeni- 
tor star mass, but depend on the local environ- 
ment. Further, the spectral behavior of the cos- 
mic rays might vary depending on the local en- 
vironment. And, finally, the average observed 
density of the galaxy might not represent the 
average density in which the secondary pho- 
tons are produced. 



3.3. A correlation between radio and 
gamma-ray emission 

The non-thermal radio emission at GHz fre- 
quencies is expected to come from synchrotron 
radiation of electrons accelerated at supernova 
shock fronts. The gamma-ray emission, on the 
other hand, is believed to have the same ori- 
gin, following the interpretation of the signal 
as the interaction of cosmic rays with the am- 
bient medium. Thus, a correlation between the 
radio and gamma-ray luminosities is expected 
to be present. Figure |2] shows the gamma-ray 
luminosity for the six objects versus their radio 
luminosity. A clear trend of increased gamma- 
ray emission at enhanced radio emission is ob- 
served. While the performed power-law fit is 
compatible with a linear correlation of the two 
emission features, the lack of statistics does 
not allow for the formulation of a quantitative 
statement on this correlation. This might still 
be a first hint that there is a connection between 
the two wavelengths, which could further be 
used to study acceleration processes. 

3.4. IHigli-energy neutrinos 

Due to the co-production of high-energy pho- 
tons and neutrinos, hadronically produced 
gamma-rays are always accompanied by a 
neutrino flux. Due to their low interaction 
probability, the detection of neutrinos requires 
the use of kilometer- scale na t ural w ater or 
ice reservoirs, see e.g. iBecker I (|2008|) . As of 
December 18, 2010, the first cubic -kilometer 
scale neutrino detector IceCube was completed 
at the geographic South Pole. The detection 
technique allows for the observation of the en- 
tire northern hemisphere at a duty cycle of 




10' 10^ 

Ughz [10" erg/s] 

Fig. 2. Luminosities at gamma-ray energies as 
measured by Fermi versus radio luminosities at 
5 GHz. Shown are the LMC, the Milky Way, M82, 
NGC253, NGC106 8 and NGC49 4 5, usi ng the data 
as summarized by iLenain et al. I ( l201lh . The line 
represents a power-law fit to the data. Due to the 
small number of sources, the error of the correlation 
is still too large to make a quantitative statement. 
Qualitatively, the result is compatible with a linear 
correlation. 



more than 99%. The main background are 
neutrinos produced in the Earth's atmosphere. 
Point source searches can be performed for 
single sources as well as by the stacking of 
a source list. If a contribution from unresolv- 
able sources is expected, which is the case 
for instance for several classes of active galac- 
tic nuclei, a search for a diffuse flux can be 
performed. Since astrophysical sources are ex- 
pected to produce relatively flat neutrino en- 
ergy spectra for interactions at the source (~ 
E~^ - E^^-^), the explicit search for enhanced 
emission at the highest energies reduces the 
background of the very steep atmospheric neu- 
trino flux (~ E^^-^). 

Concerning the case of starburst galaxies, 
the gamma-ray detection from M82 gives a 
first concrete test case of the expected neu- 
trino flux from a point source. It turns out 
that the flux of M82 as a single source is rel- 
atively low and it is not expected to be de- 
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tected within the first years of operation with 
IceCube. It may still be possible to observe 
M82 after a longer period of measurement, 
as the lifetime of IceCube is expected to be 
longer than 10 years. On the other hand, the 
stacking of a source catalog of nearby starburst 
galaxies strongly improves the detection sig- 
nificance. As the significance roughly scales 
with the signal over the square root of the back- 
ground, adding more signal to the search helps 
to reveal the signal over background. A first 
stacking search for starburst galaxies has been 
performed with partially comple ted configura- 
tions of the IceCube detector, see lDrever et al. I 
(2010); Abbasi et al. (2010). It is expected that 
within the next few years, the limits can be im- 
proved using a detector of more than twice the 
size. The detection of a diffuse flux from star- 
burst galaxies is rather challenging, since the 
maximum neutrino energy may be as low as 
10'^ eV. 



4. Molecular ions as cosmic ray 
tracers 

While protons with GeV energies and above 
contribute to the neutrino- and photon flux of 
a galaxy, low-energy cosmic rays in the keV- 
GeV range ionize the interstellar medium. In 
the Milky Way, the average ionization level 
is observe d to be of the ord e r of t!f ^ 2 ■ 
IQ-'^ s-i (iGerinetal.l 120101: iNeufeld et al.1 
I2OIO I). Hydrogen ionization immediately leads 
to the formation of Hi,, which in turn initi- 
ates the formation of larger molecules like H^, 
0H+, //2O+ etc., see e.g. (iBlacklll998l) and 
references therein. Those molecules can be ob- 
served by detectors like Herschel and ALMA 
and can be used as direct tracers of cosmic 
ray ionization. The search for molecular ions 
at potential cosmic ray acceleration sites with 
suitable targets might help to trace the sources 
of cosmic rays and by that improve the under- 
standing of the role of cosmic rays in the dy- 
namical processes of galaxies. 

In the Milky Way, several systems of su- 
pernova remnants and molecular clouds (SNR- 
MC systems) have been detected at gamma- 
ray energies in the past years. In partic- 
ular, the detections of the sources W51C 



(!Abdo et al. "2009^, W44 ("Abdo et al. "2010c'), 
W28 (A bdo et al. 2010b), IC443 (Abdo et al.] 
I2010dh and W49B ( Abdoet al.ll2010el) with 
the Fermi Gamma-ray Space Telescope are 
best-fit with a hadronic interaction model. The 
observed gamma-ray spectra can be used to es- 
timate the primary cosmic ray spectra at the 
source above GeV energies. Extrapolating the 
spectrum down to below GeV energies then 
gives the opportunity to perform calculations 
concerning interaction of the low-energy part 
of the cosmic ray spectrum leading to the ion- 
ization of the local mediu m. Details of t he cal - 
culation are presented in Becker et al. I (l201lh . 
Due to the enhanced flux of cosmic rays, the 
ionization level is expected to be enhanced by a 
few orders of magnitude at the discussed SNR- 
MC systems. In such an environment, the de- 
tection of line emission spectra from Hj and 
H3 would be the most direct tracer for cosmic 
ray ionization. 

Line emission from H^ has been dis- 
cussed previously in the context of cosmic 
ray ionization, and it represents one of the 
molecul es obser ved in astrophysical contexts, 
see e.g. dB lack I [1998) and references therein. 
With the launch of the Herschel telescope, de- 
tailed observations of the abundance of H1^ is 
now possible i n astrophysical environments. In 
llndriolo et al.l(l2010i) . for instance, an H^ abun- 
dance corresponding to an ionization rate of 
~ 2 • 10"'^/s was observed. 

Although Hj is the first product of cos- 
mic ray ionization, it is usually destroyed too 
quickly to produce significant line emission. 
However, in an environment of extreme ioniza- 
tion level as it seems to be the case for SNR- 
MC systems, the detection of Hj seems to be 
possible and would provide a unique method to 
trace the sources of cosmic rays. 

Due to their high star formation rate, star- 
burst galaxies are bound to host a larger num- 
ber of SNR-MC systems. In the future, it 
would therefore be interesting to try to com- 
bine gamma-ray measurements with the search 
for molecular line emission in order to be able 
to pin-point the cosmic ray component of the 
galaxies. 
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flux and the density of the galaxy, and not to a 
scaling with the total gas mass. 

JULIA BECKER: Assuming a constant total 
energy budget for all SNRs, this is correct, the 
scaling should be with the target density of the 
interaction. However, assuming a constant cos- 
mic ray energy density gives a scaling with the 
total gas mass. The difficulty with both argu- 
ments is that there still are a lot of assumptions 
in both statements: For instance, in both calcu- 
lations, all SNRs are assumed to have the same 
cosmic ray injection spectrum. This is not the 
case in reality. In addition, the average density 
of the galaxy does not necessarily represent the 
actual density relevant for the interactions. 

WOLFGANG KUNDT: How do you avoid 
the production of hard electrons from ioniza- 
tion? 

JULIA BECKER: The ionization cross sec- 
tions decreases rapidly with energy and can 
thus be neglected above GeV energies. This is 
why below GeV energies, ionization processes 
dominate and above GeV energies, hadronic 
interactions are prevaiUng. 



DISCUSSION 



ARNON DAR: Theoretical arguments rather 
point to a correlation between the gamma-ray 



